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ABSTRACT 

Context. Until now, there have only been seven quasars at z > 4.5 whose the high-resolution radio structure had been studied in detail 
with Very Long Baseline Interferometry (VLBI) imaging. 

Aims. We almost double the number of VLBI-imaged quasars at these high redshifts with the aim of studying their redshift-dependent 
structural and physical properties in a larger sample. 

Methods. We observed five radio quasars (J0813+3508, Jl 146+4037, J1242+5422, J161 1+0844, and J1659+2101) at 4.5 < z < 5 
with the European VLBI Network (EVN) at 1.6 GHz on 29 October 2008 and at 5 GHz on 22 October 2008. The angular resolution 
achieved ranges from 1.5 to 25 milli-arcseconds (mas), depending on the observing frequency, the position angle in the sky, and the 
source's celestial position. 

Results. The sources are all somewhat extended on mas scales, but compact enough to be detected at both frequencies. With one 
exception of a flat-spectrum source (J161 1+0844), their compact emission is characterised by a steep radio spectrum. We found no 
evidence of Doppler-boosted radio emission in the quasars in our sample. The radio structure of one of them (J08 1 3+3508) is extended 
to ~7", which corresponds to 43 kpc projected linear size. Many of the highest redshift compact radio sources are likely to be young, 
evolving objects, far-away cousins of the powerful gigahertz peaked-spectrum (GPS) and compact steep-spectrum (CSS) sources that 
populate the Universe at lower redshifts. 

Key words, radio continuum: galaxies - galaxies: active - quasars: general - techniques: interferometric 

1 . Introduction (i) their radio flux density is relatively low for VLBI, requiring 

high data rate phase-referenced observations; (ii) due to their 

Quasars at the highest redshifts place strong constraints on the steep spectra, they are generally not considered useful for VLBI 

early cosmological evolution of active galactic nuclei (AGNs) experiments; (Hi) for most of them, spectroscopic redshifts are 

and the growth of their central supermassive (~10 9 M ) black simply not yet available. Indeed, the fine-scale radio structure of 

holes. The AGN activity obser ved at h igh z indicates that feed- J0836+0054 (z = 5.77) and J1427+3312 (z = 6.12) could fortu- 

back processes (e.g. Best et al. [2005 ) may have played an im- nately be revealed because their record-breaking redshifts made 

portant role in the early galaxy and cluster evolution. The ulti- these sources attractive for high-resolution VLBI imaging. 

Ci ' mate evidence for AGN iets is provided by Very Long Baseline ™ , 7T nT ■ c ,. , . , ,. c ^ , ,. , , 

i T „ ■} " • , ,• , • The VLBI images of the highest-redshift, known radio-loud 

Interferometry (VLBI) observations in the radio domain. T1 .~ n . , ■ „f T .. „ , onn< n , 

J v ' quasar J1427+3312 at 1.4 GHz (Momjian et al. 20Q8J and 

Compact radio sources have flat spectra (with power-law 1.6 GHz (Frey et al. 2008b) reveal a prominent double structure. 

spectral index a > -0.5; S cc v", where S is the flux den- The two resolved components separated by ~28 mas (~160 pc) 

sity and v the frequency) from the synchrotron self-absorption resemble a compact symmetric object (CSO, Wilkinson et al. 

(Kellermann & Paulinv-Toth I1969i l. In the past, VLBI targets 1994). (To calculate linear sizes and luminosities, we assume a 

were traditionally selected by their flat overall radio spectrum to flat cosmological model with Ha = 70 km s _1 Mpc~' , Q m = 0.3, 

ensure detectability. Recent surveys that did not apply spectral and Qa = 0.7 throughout this paper.) CSOs are a class of very 

selection criteria (e.g. Mosoni et al. 120061 Frev et al. l20~08al) . and young (< 10 4 yr) sources typically found in radio galaxies at 

VLBI observations of individual z~6 quasars (J0836+0054, Frey much lower redshifts (z < 1). Apart from the structural similar- 

et al. 120031120051 J1427+3312,Frey etal.|2008b; Momjian et al. ity to CSOs, there are several indications of the youthfulness of 

2008) indicate that there is a less-known steep-spectrum popu- this quasar: the steep radio spectrum coupled with the compact 

lation of compact radio AGNs at high redshifts. These sources structure, the broad absorption lines, and the possible intrinsic 

have so far been able to escape discovery for a variety of reasons: X-ray absorption. 
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Table 1. Quasars at z > 4.5 imaged earlier with VLBI by in- Table 2. Our VLBI targets, five SDSS/FIRST radio quasars at 
creasing redshift. z > 4.5, and a close radio companion of one of them. 



Source name z v Network Peak Ref. 



GHz mJy/beam 



J1235-0003 


4.69 


1.4 


VLBA 


17 


1 


J1205-0742N 


4.70 


1.4 


VLBA 


0.2 


2 


J1205-0742S 


4.70 


1.4 


VLBA 


0.2 


2 


J 1430+4204 


4.72 


5 


EVN 


161 


3 






15 


VLBA 


159 & 200 


4 


J1451-1512 


4.76 


5 


EVN 


50 


5 


J0913+5919 


5.11 


1.4 


VLBA 


19 


1 


J0906+6930 


5.47 


15 


VLBA 


115 


6 






43 


VLBA 


42 


6 


J0836+0054 


5.77 


1.6 


EVN 


0.8 


7 






5 


EVN 


0.3 


8 


J1427+3312 


6.12 


1.4 


VLBA 


1.0 


9 






1.6 


EVN 


0.5 


10 






5 


EVN 


0.2 


10 



Notes: Col. 1 - source name (J2000); Col. 2 - spectroscopic redshift; 
Col. 3 - observing frequency (GHz); Col. 4 - interferometer array 
(EVN: European VLBI Network, VLBA: Very Long Baseline Array); 
Col. 5 - peak brightness (mJy/beam); Col. 6 - references (1: Momjian 
et al. [2004 2: Momjian et 31. 120051 3: Paragi et al. [T999l 4: Veres et 
al. 120101 5: L.I. Gurvits et al., in preparation; 6: Romani et al. 120041 7: 
Fr ey et al.l2003l 8: Frey et al. l20031 9: Momjian et al. 2008} 10: Frey et 
al. l2008bt . 



To our knowledge, there are only seven radio-loud AGNs 
at z > 4.5 that have been imaged with VLBI prior to our ex- 
periment reported here (Table [TJ. Four of them (J0906+6930, 
J 1 235 -0003, J 1430+4204, and J 145 1 - 1 5 1 2) are compact, prac- 
tically unresolved flat-spectrum radio sources, while others 
(J0836+0054, J0913+5919, and J1427+3312 with a double 
structure) have a compact or somewhat resolved appearance and 
a steep radio spectrum in the GHz frequency range (in the ob- 
server's frame). The double quasar J1205-0742 also listed in 
the table is a special case from our point of view, with mea- 
sured brightness temperatures indicating extreme nuclear star- 
bursts rather than radio-loud AGNs (Momjian et al. 2005). 

Our goal was to substantially increase the number of radio- 
loud AGNs at z > 4.5 imaged with VLBI. One could expect 
to find either more "classical" core-jet sources or other steep- 
spectrum quasars, possibly with CSO-like double structures. If 
the latter are found, in a decade-long term, VLBI monitoring 
would eventually allow measurements of the component expan- 
sions and thus facilitate direct estimations of the kinematic age 
of the sources. Otherwise, the compact core-jet sources are po- 
tentially valuable additions for comparing of the mas-scale struc- 
tures of quasars at low and high redshift. Efforts to use classi- 
cal cosmological tests - the angular size-redshift relation (e.g. 
Gurvits et al. 1999) and the apparent proper motion-redshift re- 
lation (e.g. Kellermann et al. 1 19991 1 - would also benefit from 
data on an increased sample of extremely distant quasars, since 
the predictions of the various cosmological world models are 
markedly different at the highest redshifts. Despite the prac- 
tical difficulties, there is continuous interest in these tests in 
the community (e.g. Sahni & Starobinsky 2000; Vishwakarma 
120001 1200T1 Lima & Alcaniz[2002 Chen & Ratra|2003 Jackson 
[2004ll2008T >. 

In this paper, we report on our dual-frequency VLBI imaging 
observations of five previously unexplored radio-loud quasars at 
4.5 < z < 5. Our sample selection method is described in Sect. [2] 
The experiments and the data reduction are explained in Sect. [3] 
The observed radio properties of the sources are given in Sect.|4] 



Source coordinates 


z 


r 


FIRST peak 


S 1.4 








mJy/beam 


mJy 


08 13 33.32 +35 08 10.8 


4.92 


20.8 


23.2 


25.2 


08 13 32.89 +35 08 14.9* 






11.8 


11.8 


11 46 57.79 +40 37 08.6 


5.01 


21.0 


12.5 


12.5 


12 42 30.58 +54 22 57.3 


4.73 


20.9 


19.7 


20.2 


16 11 05.64 +08 44 35.4 


4.54 


19.7 


8.8 


8.8 


16 59 13.23 +21 01 15.8 


4.78 


21.5 


28.7 


28.8 


Notes: Col. 1 - a priori source J2000 right ascension ( h m s ) and declina- 



tion (° ' ") from SDSS; Col. 2 - spectroscopic redshift; Col. 3 - SDSS 
r magnitude (Schneider et al. I2007t ; Col. 4 - FIRST peak brightness 
(mJy/beam); Col. 5 - FIRST integral 1.4-GHz flux density (mJy). 
* optically unidentified radio companion of the previous source; coordi- 
nates from FIRST 

Our results are discussed in Sect. [5] Conclusions are drawn in 
Sect. [6] 

2. Target selection 

For the high-resolution VLBI observations, we choose five 
quasars from the Sloan Digital Sky Survey (SDSS) Data Release 
5 (DR5) quasar catalogue (Schneider et al. 2007). All of them 
are identified with an unresolved (< 5") radio source in the Very 
Large Array (VLA) Faint Images of the Radio Sky at Twenty- 
centimeter (FIRST) survejQ (White et al. [19971 ), with 1.4-GHz 
total flux densities 8.8 mJy < S 1.4 < 28.8 mJy. These five 
quasars, together with J0913+5919 (already studied with VLBI 
by Momjian et al. 2004; see Table [TJ, are the only such radio 
quasars in the northern hemisphere in the Schneider et al. ( 2007 ) 
catalogue with z > 4.5 and S 1.4 > 5 mJy. Important is that the 
(otherwise unknown) radio spectral index of the sources was not 
used as a selection criterion. 

One of the sources (J0813+3508, z = 4.92) has a close 
(< 7"), optically unidentified radio companion in the FIRST 
catalogue, which appears compact on an arcsecond scale in the 
1.4-GHz VLA FIRST image (Fig.Q). We also included this ob- 
ject as the sixth target in our VLBI experiment, in the hope that 
any possible relation between the two apparently nearby sources 
could be studied. The basic parameters of our target sources are 
listed in Table [2] 

Our recent experience with a larger sample of SDSS/FIRST 
quasars (Mosoni et al. 2006; Frey et al. 2008a) indicates that 
the sources identified as optical quasars and unresolved FIRST 
objects with S 1.4 > 20 mJy have a nearly 90% chance to be suc- 
cessfully detected with the European VLBI Network (EVN) at 
5 GHz, using the SDSS coordinates as a priori values. Although 
most of the sources are weaker in the present sample, we used 
longer integration times in the VLBI experiment to ensure safe 
detections. 

3. VLBI observations and data reduction 

We observed our targets with the EVN at 1 .6 GHz on 29 October 
2008 and at 5 GHz on 22 October 2008. Due to the high red- 
shifts, the observed frequencies correspond to ~9— 10 GHz and 
~28-30 GHz in the rest frame of the quasars. At a recording 
rate of 1024 Mbit s _1 , ten antennas of the EVN participated 
in the experiment at 1.6 GHz: Effelsberg (Germany), Jodrell 

1 http://sundog.stsci.edu 
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08 13 35.0 34.5 34.0 33.5 33.0 32.5 32.0 

RIGHT ASCENSION (J2000) 

Fig.l. The 1.4-GHz image o f J081 3+3508 from the VLA 
FIRST survey (White et al. 1997). The peak brightness 
is 22.5 mJy/beam, the lowest contour levels are drawn at 
+0.25 mJy/beam, the positive contour levels increase by a factor 
of 2. The circular Gaussian restoring beam size is 5'.'4 (FWHM), 
as indicated in the lower-left corner. The northwestern radio 
companion (top right) is catalogued as another source in FIRST, 
6'.'8 away from the quasar (in the centre). 



Bank Mk2 (UK), Medicina, Noto (Italy), Toruh (Poland), Onsala 
(Sweden), Sheshan, Nanshan (RR. China), Robledo (Spain), and 
the phased array of the 14-element Westerbork Synthesis Radio 
Telescope (WSRT, The Netherlands). The eight participating an- 
tennas in the 5-GHz experiment were Effelsberg, Jodrell Bank 
Mk2, Medicina, Noto, Toruh, Sheshan, Nanshan, and the WSRT. 
On both days, the observations lasted for a total of 1 1 h. Eight 
intermediate frequency channels (IFs) were used in both left and 
right circular polarisations. The total bandwidth was 128 MHz 
per polarisation. The correlation of the recorded VLBI data took 
place at the EVN Data Processor at the Joint Institute for VLBI 
in Europe (JIVE), Dwingeloo, the Netherlands. 

All the target sources were observed in phase-reference 
mode. This allows us to increase the coherent integration time 
spent on the source and thus to improve the sensitivity of the 
observations. Phase-referencing involves regularly interleaving 
observations between the target source and a nearby bright and 
compact reference source (e.g. Beasley & Conway 1995). The 
delay, delay rate, and phase solutions derived for the phase- 
reference calibrator were interpolated and applied for the re- 
spective target within the cycle time of ~7 minutes. The tar- 
get sources were observed for ~4.5-minute intervals in each 
cycle. The exception was the pair of sources, J0813+3508 and 
its optically unidentified apparent companion (Table |2l Fig.[T|). 
They were observed using the same phase-reference source, for 
~2.5 minutes subsequently in each cycle. To achieve nearly the 
same total on-source integration time as for the other targets 
(~60 min), more cycles were scheduled for this pair. 

The suitable phase-reference calibrator sources 
(J0815+3635, Jl 146+3958, J1253+5301, J1608+1029, and 
J1656+1826) were selected from the Very Long Baseline Array 
(VLBA) Calibrator SurveyQ. The angular separations between 
the calibrators and the corresponding targets range from 0°64 



http : //www . vlba . nrao . edu/astro/calib/ index . shtml 



to 2°65 (Table [3]). The positional uncertainties of the calibrators 
in the International Celestial Reference Frame (ICRF) are 
0.4-4.3 mas. 

The US National Radio Astronomy Observatory (NRAO) 
Astronomical Image Processing System (AIPS; Diamond 1995) 
was used for the data calibration. The visibility amplitudes 
were calibrated using the antenna gains, and the system tem- 
peratures regularly measured at the antennas during the ex- 
periments. Fringe-fitting was performed for the five calibrators 
mentioned above, and the fringe-finder sources (J0555+3948, 
J0927+3902, Jl 159+2914, and J1331+3030) using 3-min solu- 
tion intervals. The data were exported to the Caltech Difmap 
package (Shepherd et al. 119941) for imaging. The conventional 
mapping procedure involving several iterations of CLEANing 
and phase (then amplitude) self-calibration resulted in the im- 
ages and brightness distribution models for the calibrators. 
Overall antenna gain correction factors (~10% or less) were de- 
termined and applied to the visibility amplitudes in AIPS. Then 
fringe-fitting was repeated in AIPS, now taking the clean com- 
ponent models of the phase-reference calibrator sources into ac- 
count. The residual phase corrections resulted from their non- 
pointlike structure were considered this way. The solutions ob- 
tained were interpolated and applied to the target source data. 
The visibility data of the target sources, unaveraged in time and 
frequency, were also exported to Difmap for imaging. The nat- 
urally weighted images at 1 .6 GHz and 5 GHz (Fig. |2]i were 
made after several cycles of CLEANing in Difmap. Phase-only 
self-calibration was applied for the brighter sources (i.e. when 
the sum of the CLEAN component flux densities exceeded 
~10 mJy) over time intervals not shorter than the length of scans 
spent on the sources (see Table [3). In the cases where phase 
self-calibration was not attempted at all, we expect a loss of co- 
herence of about 5% (cf. Marti- Vidal et al. 2010) in the phase- 
referencing process, which may lead to an underestimate of the 
flux density values by this factor. In the images, the lowest con- 
tours are drawn at ~3<x image noise levels. The expected theoret- 
ical thermal noise values were 15-22 /iJy/beam (lcr), assuming 
no data loss during the experiment. The image parameters are 
summarised in Table [3] 

4. Results 

The five targeted high-redshift quasars were clearly detected and 
imaged with the EVN at both 1.6 GHz and 5 GHz (Fig. 0. 
However, we could not detect any mas-scale compact radio 
source at the position of the nearby (6'.'8) FIRST radio compan- 
ion (Fig.[T} to J0813+3508, at a brightness level of 0.3 mJy/beam 
(3<x). This indicates that the secondary source is extended and 
thus completely resolved in our VLBI experiments; notably, the 
position angle of the jet-like structure in our VLBI images of 
J0813+3508, especially at 1.6 GHz (Fig. |2] top left) practically 
coincides with that of the FIRST companion. Based on this di- 
rectional coincidence, we believe that there is a physical asso- 
ciation between the two radio sources. This is further supported 
by the fact that the visibility amplitudes on the shortest and most 
sensitive baseline from Effelsberg to Westerbork were somewhat 
higher than could be fully accounted for with the CLEAN model 
components derived from the full VLBI array data. From poorly 
constrained Gaussian model fitting, there was a hint on an ex- 
tended component about 076 away from the compact "core" to 
the northwest in the continuation of the ~10-mas scale jet and in 
the direction of the FIRST radio companion. 

The coordinates of the brightness peaks are estimated from 
the 5-GHz images using the AIPS task maxfit and listed in 
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Table 3. VLBI image parameters for Figj2] 



Source name cj> v Peak Contour Restoring beam 
° GHz mJy/beam masxmas PA (°) 



J0813+3508 


1.50 


1.6 


11.4 


0.20 


19.0x4.9 


5.3 






5 


4.0* 


0.15 


5.4 x 1.5 


8.6 


J 1146+4037 


0.64 


1.6 


14.7 


0.30 


20.1 x4.4 


5.4 






5 


6.0* 


0.20 


5.1 x 1.5 


3.6 


J1242+5422 


2.09 


1.6 


12.1 


0.60 


17.1 x4.5 


22.2 






5 


7.9 


0.15 


4.5 x 1.5 


22.7 


J161 1+0844 


1.83 


1.6 


11.1 


1.00 


25.7 x 3.4 


9.6 






5 


9.4 


0.25 


7.1 x 1.5 


9.1 


J1659+2101 


2.65 


1.6 


16.6 


0.30 


22.4 x 3.7 


11.6 






5 


4.2 


0.50 


5.8 x 1.4 


11.7 



Notes: Col. 1 - source name (J2000); Col. 2. - angular separation from the phase-reference calibrator source (°); Col. 3 - observing frequency 

(GHz); Col. 4 - peak brightness (mJy/beam); cases where phase self-calibration was not applied are marked with asterisks; Col. 5 - lowest 
contour level (mJy/beam) corresponding to ~3cr image noise; Col. 6 - Gaussian restoring beam size (masxmas); Col. 7 - restoring beam major 

axis position angle (°) measured from north through east. 



Table 4. Measured coordinates of the five radio quasars at z > 
4.5 imaged with VLBI. 



Source name 


Right ascension 


Declination 


Uncertainties 




h m s 


O 1 If 


mas 


mas 


J0813+3508 


08 13 33.32789 


35 08 10.7698 


0.4 


0.5 


Jl 146+4037 


11 46 57.79043 


40 37 08.6256 


0.3 


0.4 


J1242+5422 


12 42 30.58994 


54 22 57.4524 


3.9 


1.9 


J161 1+0844 


16 11 05.65000 


08 44 35.4776 


0.5 


0.4 


J1659+2101 


16 59 13.22857 


21 01 15.8087 


1.0 


1.1 



Notes: Col. 1 - source name; Col. 2 - J2000 right ascension ( h m s ); 
Col. 3 - J2000 declination (° ' "); Col. 4 & 5 - estimated positional 
uncertainties in right ascension and declination (mas). 



Table[4] The uncertainties are determined by the phase-reference 
calibrator source position accuracy, the target-calibrator angular 
separation, the angular resolution of the interferometer array, and 
the signal-to-noise ratio. 

Difmap was used to fit circular Gaussian brightness distribu- 
tion model components to the interferometric visibility data of 
the five quasars detected. We calculated the two-point spectral 
indices for the bright central components from the model flux 
densities. From the 5-GHz data, we also derived the brightness 
temperatures of these dominant components as 

r B = 1.22xl0 12 (l+z)^-2 [K], (1) 

where z is the redshift, S the flux density (Jy), v the observ- 
ing frequency (GHz), and i? the full width at half maximum 
(FWHM) size of the Gaussian measured in mas (e.g. Condon 
et al. 1 1 9821 ). The model parameters, and the values derived from 
them are given in Table [5] 

5. Discussion 

All of the z > 4.5 quasars we imaged are somewhat resolved, of- 
ten with structures extending up to several 10 mas (Fig. [2]). For 
comparison, a 1 mas angular size corresponds to 6.3-6.6 pc pro- 
jected linear size, depending on the actual redshift of the object. 
The measured brightness temperatures (rB~10 8 -10 9 K; Table|5j 
clearly indicate AGN activity, but are at least an order of magni- 
tude lower than the equipartition value estimated for relativistic 
compact jets (7/ B , e q - 5 X 10 10 K; Readhead fl994| l. It may sug- 
gest that (0 in the radio spectrum of the sources we probe regions 



far away from the peak frequency caused by synchrotron self- 
absorption; (zz) the jet viewing angles are at least moderate, and 
in fact we experience Doppler-deboosting of the emission; or 
(Hi) the intrinsic brightness temperatures of compact jets at very 
high redshift are significantly lower than the equipartition value. 
Except for one flat-spectrum source (J 161 1+0844), the spectra 
of the quasar "cores" (i.e. the main components) in our sam- 
ple are steep (a < -0.5; Table|5]). The intrinsic brightness tem- 
peratue was determined for the blazar J 1430+4204 (z = 4.72) 
by comparing the variability brightness temperature and the one 
measured with VLBI. It was found that T^int - ?B,eq (Veres et 
al. l2010l l. so there is no reason to assume that the intrinsic bright- 
ness temperatures are generally lower at z > 4.5. 

A possible signature of the blazar nature of a source is its 
strong variability. We can look for an indication of long-term 
variability, if we compare the 1 .4-GHz flux densities taken from 
the VLA FIRST catalogue (Table |2 last column) with our 1.6- 
GHz flux densities measured in our EVN run (Table|5]). Since our 
angular resolution is much higher, if a source is stationary, then 
one would expect the VLBI flux densities to be lower than or at 
best equal to the VLA values within the uncertainties. This is the 
case for all our sources but J 1 6 1 1 +0844. That J 1 6 1 1 +0844 - the 
only flat-spectrum source in our sample - has nearly 50% higher 
flux density measured at a later epoch at much higher resolu- 
tion proves its non-stationarity. Both the variability and the flat 
spectrum argue for Doppler boosting, but the measured bright- 
ness temperature of this source does not support it. In the case of 
J161 1+0844, the observed one-sided compact core-jet structure 
at 5 GHz and the quasar classification can just be made con- 
sistent with the simple orientation-dependent unified picture of 
radio-loud AGNs (e.g. Urry & Padovani 1995). The apparent de- 
boosting (Doppler factor 5 = Tb/Tb,^ - 0.1) can be explained 
if, e.g., we assume a = 45° jet angle to the line of sight and an 
extremely (but not impossibly) large bulk Lorentz factor (y^35). 

The case of J0813+3508 is similarly puzzling. If, as our data 
suggest, its FIRST companion (Fig. [TJ) is associated with the 
quasar at z = 4.92, the projected linear size of the source is 
43 kpc. This means that J0813+3508 could be the quasar with 
the most extended radio jet known at an extremely high red- 
shift (cf. Cheung et al. IIOTBI 120081 . The markedly one-sided 
arcsecond-scale radio structure suggests Doppler boosting and 
a small jet angle to the line of sight. However, the measured 
~10 9 K brightness temperature and the steep spectrum of the 
compact VLBI component do not support this view. It is pos- 
sible that there is a significant misalignment between the mas- 
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J0813+3508 




1.6 GHz 



J0813+3508 




5 GHz 



Right Ascension (mas) 
~l 1 1 1 1 i 1 1 

J1146+4037 



10 -10 

Right Ascension (mas) 



1.6 GHz 



J1146+4037 




5 GHz 



Right Ascension (mas) 



-10 
Right Ascension (mas) 



J1242+5422 




1.6 GHz 



J1242+5422 




5 GHz 



Right Ascension (mas) 





Right Ascension (mas) 



Fig. 2. The naturally weighted 1.6-GHz (left column) and 5 -GHz (right column) VLBI images of the quasars. The image parameters 
(peak brightness, lowest contour level corresponding to ~3cr image noise, restoring beam size, and orientation) are listed in Table[3] 
The full width at half maximum (FWHM) of the Gaussian restoring beam is indicated with an ellipse in the lower-left corners. The 
positive contour levels increase by a factor of 2. The coordinates are related to the brightness peak of which the absolute position is 
given in Table [4] 
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Fig. 2. continued. 



Table 5. The parameters derived for the quasars from VLBI imaging. 



Source name 


Complete VLBI flux density (mJy) 


Main component flux density (mJy) 


a 


Size 


T B 


L 




at 1.6 GHz 


at 5 GHz 


at 1.6 GHz 


at 5 GHz 




mas 


10 9 K 


10 26 W Hz" 1 


J0813+3508 


17.1+0.9 


7.3+0.4* 


12.8+0.6 


6.8+0.3* 


-0.6 


1.18+0.02 


1.5+0.1 


8.4+0.4 


Jl 146+4037 


15.5±0.8 


8.6+0.4* 


15.5+0.8 


8.4+0.4* 


-0.5 


0.74+0.01 


4.5+0.3 


10.7+0.5 


J1242+5422 


17.7+0.9 


9.7+0.5 


17.1+0.9 


9.4+0.5 


-0.5 


0.67+0.01 


5.9+0.5 


9.0+0.5 


J161 1+0844 


13.0+0.8 


12.9+0.6 


13.0+0.8 


12.4+0.6 


-0.0 


0.85+0.01 


4.7+0.3 


4.6+0.2 


J1659+2101 


29.3+1.5 


10.6+0.7 


19.7+1.0 


10.6+0.7 


-0.6 


3.07+0.12 


0.3+0.04 


12.3+0.8 



Notes: Col. 1 - source name; Col. 2-3 - complete VLBI flux density at 1.6 GHz and 5 GHz (mJy); Col. 4-5 - main component VLBI flux density 
at 1.6 GHz and 5 GHz (mJy); Col. 6 - two-point spectral index of the main component; Col. 7 - fitted circular Gaussian diameter (FWHM) for 
the main component at 5 GHz (mas); Col. 8 - rest-frame brightness temperature at 5 GHz (10 9 K); Col. 9 - rest-frame monochromatic 5-GHz 
luminosity (10 26 W Hz~'). The VLBI flux density calibration uncertainties are assumed as 5%. The statistical errors of the fitted model parameters 
are estimated according to Fomalont d 1999b . 

* flux density value corrected for an estimated coherence loss of 5% 



scale and arcsecond-scale jet, or - less probably - the jet in this 
quasar is intrinsically one-sided. An alternative scenario is that 
we see an expanding and extended radio source with double hot 
spots, of which the nearest and approaching one is detected with 
VLBI. Future sensitive radio interferometric observations at in- 



termediate angular resolutions could help reveal the true nature 
of this object (or these objects). 

The remaining three sources (Jl 146+4037, J1242+5422 and 
J1659+2101) are quite alike in their compact structure (Fig. 13, 
steep spectra, and luminosity (Table|5]). The quasar J1659+2101 
is somewhat more resolved than the others. Notably, the two 
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most distant (z~6) quasa rs (J083 6+0054, Frey et al. [20 031 [2005 1 
J1427+3312, Frey et al. 2008b Momjian et al. 120081 observed 
with VLBI to date share much similar properties. 

The rest-frame 5-GHz luminosities of the main components 
of our sources (L~10 26 - 10 27 W Hz -1 ; Tabled are compara- 
ble to the values of a large sample of gigahertz peaked-spectrum 
(GPS) and compact steep-spectrum (CSS) sources compiled by 
O'Dea 019981 1. Our values refer to the VLBI components and can 
therefore be regarded as lower limits to the radio luminosities of 
the whole sources. We found no evidence of any significantly 
Doppler-boosted radio emission in our cases, so that the mea- 
sured luminosities indicate intrinsically powerful sources. 

GPS and CSS sources are thought to be young, evolving ob- 
jects, and/or perhaps "frustrated" ones that are confined by the 
dense ambient gas (see O'Dea 1998 for a review). A subclass 
of these sources are the CSOs, which show nearly symmetric 
double or triple morphology when imaged at VLBI resoulution. 
According to the model of Falcke et al. (2004), our high-redshift 
steep-spectrum objects may represent GPS sources at early cos- 
mological epochs. The first generation of supermassive black 
holes could have powerful jets that developed hot spots well in- 
side their forming host galaxy, on linear scales of 0.1-10 kpc. 
Taking the relation between the source size and the turnover fre- 
quency observed for GPS sources into account, and for hypothet- 
ical sources matching the luminosity and spectral index of ours, 
Falcke et al. ([2004) predict that the angular size of the smallest 
(~ 100 pc) of these early radio-jet objects would be in the order of 
10 mas at z~5, and the observed turnover frequency in their ra- 
dio spectra would be around 500 MHz. These angular sizes are 
indeed seen in our VLBI images (Fig. To confirm the low- 
frequency turnover, high-resolution interferometric flux density 
measurements would be needed at multiple frequencies down to 
the 100-MHz range - a task well-suited to the Square Kilometre 
Array (SKA) now under development, in its high-resolution con- 
figuration. In a sense, the observed spectral properties of high- 
redshift sources are a resurrection of the indication on "humped" 
spectra of high-redshift quasars found more than 20 years ago 
when the highest known redshift was below 4 (Peterson et al. 
[19821 O'Dea[1990ll. 

Currently we do not see many blazars at z~5 or higher. 
According to Table [1] the most distant blazar imaged with VLBI 
to date (J0906+6930; Romani et al. 120041 has a redshift of 
5.47. On the other hand, one would naively expect highly rela- 
tivistically beamed, flat-spectrum radio sources to dominate the 
high-redshift VLBI samples. Whether the apparent absence of 
very high-redshift blazars is a result of some selection effect 
(i.e. the lack of redshift measurements), the poor statistics due 
to the small sample, or has a real physical cause, is to be ad- 
dressed with further observations. Indeed, there are indications 
that blazars are not very common at the highest redshifts. Based 
on the evolving gamma-ray luminosity function, Inoue et al. 
(120101 1 expect that the Fermi Gamma-ray Space Telescope will 
find a few (i.e. the order of unity) blazars at z~6 over a period 
of 5 years. For the strong sources, there is a positive correlation 
between the gamma-ray and the parsec-scale synchrotron radio 
emission, with the gamma-detected sources having on average 
higher brightness tempeartures in the radio (e.g. Kovalev et al. 
2009). Thus the Fermi Large Area Telescope (LAT) potentially 
selects the brightest radio sources that would be natural targets 
for follow-up with VLBI. 



6. Conclusions 

We imaged five distant radio quasars at 4.5 < z < 5 with the 
EVN at two frequencies (1.6 and 5 GHz), almost doubling the 
currently known sample of quasars imaged with VLBI at z > 4.5. 
The phase-referenced observations allowed us to derive accu- 
rate astrometric positions for our targets. The slightly resolved 
mas- and 10-mas-scale radio structures, the measured moderate 
brightness temperatures (~10 8 -10 9 K), and the steep spectra in 
all but one case suggest that our sample of compact radio sources 
at z > 4.5 is dominated by objects that do not resemble blazars. 
One of the quasars (J0813+3508) is likely to be extended to ~7" 
which corresponds to a 43 kpc projected linear size. It is possible 
that we see young, evolving, and compact GPS-like objects that 
are signatures of early galaxy formation where the expanding 
powerful synchrotron radio sources are confined by the dense 
interstellar medium. 
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